COMBINATORICS OF CREMONA MONOMIAL MAPS 



ARON SIMIS AND RAFAEL H. VILLARREAL 



Abstract. One studies Cremona monomial maps by combinatorial means. Among the results 
is a simple integer matrix theoretic proof that the inverse of a Cremona monomial map is 
also defined by monomials of fixed degree, and moreover, the set of monomials defining the 
inverse can be obtained explicitly in terms of the initial data. A neat consequence is drawn 
for the plane Cremona monomial group, in particular the known result saying that a plane 
Cremona (monomial) map and its inverse have the same degree. Included is a discussion about 
the computational side and/or implementation of the combinatorial invariants stemming from 
these questions. 



1. Introduction 

The expression "birational combinatorics" has been introduced in [10] to mean the combina- 
torial theory of rational maps P n_1 p m_1 defined by monomials, along with natural integer 
arithmetic criteria for such maps to be birational onto their image varieties. As claimed there, 
both the theory and the criteria were intended to be a simple transcription of the initial geo- 
metric data. Yet another goal is to write characteristic-free results. Thus, here too one works 
over an arbitrary field in order that the theory be essentially independent of the nature of the 
field of coefficients, specially when dealing with squarefree monomials. 

In this paper, we stick to the case where m = n and deal with Cremona maps. An important 
step has been silently taken for granted in the background of [101 Section 5.1.2], namely, that 
the inverse of a Cremona monomial map is also defined by monomials. To be fair this result 
can be obtained via the method of [10| Section 3] together with the criterion of [8]; however the 
latter gives no hint on how to derive explicit data from the given ones. 

Here we add a few steps to the theory, by setting up a direct way to convert geometric 
results into numeric or combinatorial data regardless of the ground field nature. The conversion 
allows for an incursion into some of the details of the theory of plane Cremona maps defined 
by monomials. In particular, it is shown that the group of such maps under composition is 
completely understood without recurring to the known results about general plane Cremona 
maps. Thus, one shows that this group is generated by two basic monomial quadratic maps, up 
to reordering of variables in the source and the target. The result is not a trivial consequence 
of Noether's theorem since the latter requires composing with projective transformations, which 
is out of the picture here. Moreover, the known proofs of Noether's theorem (see, e.g., pQ) 
reduce to various especial situations, passing through the celebrated de Jonquieres maps which 
are rarely monomial. 

The well-known result that a plane Cremona map and its inverse have the same degree is 
shown here for such monomial maps by an easy numerical counting. The argument for general 
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plane Cremona maps is not difficult but requires quite a bit of geometric insight and preparation 
(see, e.g., HJ Proposition 2.1.12]). 

Monomial Cremona maps have been dealt with in [5] and in [6], but the methods and some 
of the goals are different and have not been drawn upon here. 

2. Tools of integer linear algebra 

Recall that if a = (ai, . . . , a n ) G M n , its support is defined as supp(a) = {i \ 04 7^ 0}. Note that 
we can write a = a + — a~ , where a + and a~ are two non-negative vectors with disjoint support. 
The vectors a + and a~ are called the positive and negative part of a respectively. Following a 
familiar notation we write \a\ = a\ + ■ ■ ■ + a n . 

The following result is essentially embodied in the treatment given in the preliminaries of [10J. 
For easy reference we chose to isolate it along with its complete proof. 

Lemma 2.1. Let vi,...,v n be a set of vectors in N n such that \vi\ = d > 1 for all i and 
det(A) = ±d, where A is the nxn matrix with column vectors v±, . . . , v n . Then A~ 1 (ei — ej) G 27* 
for all i,j. 

Proof. Fixing indices i, j, there are Ai, . . . , A n in Q such that A~ 1 {ei — ej) = Ysk=i ^fc e fc- Notice 
that A~ 1 (e i ) is the ith column of A' 1 . Set 1 = (1, ... , 1). Since 1A = dl, we get l/d = lA^ 1 . 
Therefore |A~ 1 (e i )| = |A _1 ( e i)l = l / d and Efc X k = 0. Then we can write 

n n 

A^ 1 (e i - ej) = ^2 ^k{ei - e x ) ^ - ej = ^ X k (v k - v{). 

k=2 k=2 

Thus there is 7^ s G N such that s(ei-ej) belong to ^{vi—Vk} k=2 , the subgroup of 7L n generated 
by {v\ — Vk} k=2 - By [lQl Lemma 2.2 and Theorem 2.6], the quotient group Z n j7L{v\ — Vk} k=2 is 
free, in particular has no nonzero torsion elements. Then we can write 

ei - ej = r] 2 (v 2 — v±) H h 7] n (v n - v\), 

for some 77^ 's in Z. Since 7L{v\ — Vk} k=2 is also free (of rank n — 1), the vectors v 2 — vi, . . . ,v n — vi 
are linearly independent. Thus At = % 6 Z for all k > 2, hence ultimately A~ 1 (ei — ej) G Z n . □ 

Next we state our main result of integer linear algebra nature. Its geometric translation and 
applications will be given in Section [3l 

Theorem 2.2. Let vi,...,v n be a set of vectors in N n such that \vi\ = d > 1 for all i and 
det(.A) = ±d ; where A is the nxn matrix with column vectors v\, . . . ,v n . Then there are unique 
vectors {3%, . . . ,/3 n >7 G N n such that the following two conditions hold: 

(a) A0i = 7 + &i for all i, where /3j,7 and ei are regarded as column vectors; 

(b) The matrix B whose columns are Pi, . . . ,(3 n has at least one zero entry in every row. 

Moreover, det(B) = ±(\j\ + l)/d = ±|A| for all i. 

Proof. First we show the uniqueness. Assume that P[, . . . , f3' n , 7' is a set of vectors in N n such 
that: (a') A(3[ = 7' + e% for all i, and (b') The matrix B' whose column vectors are /3[, . . . , (3' n 
has at least one zero entry in every row. Let A = (Aj) and A' = (A^) be non-negative vectors 
such that j4 _1 (7 — 7') = A' — A. Then from (a) and (a') we get 

(2.1) Pi - /3i = A- 1 ( 7 - 7 ') = A' - A, Vi =}► p ik - P' ih = A' k - A fcj Vi, k, 
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where = (fin, . . . , (3 in ) and /?• = (J3' a , /?■„). It suffices to show that A = A'. If A' k > A k for 
some k, then, by Eq. (|2.ip . we obtain 0^. > for i = 1, ... ,n, which contradicts (b). Similarly 
if A' k < Afc for some k, then, by Eq. (|2.ip . we obtain > for i = 1, . . . , n, which contradicts 
(b'). Thus A fc = A' k for all k, i.e., A = A'. 

Next we prove the existence of P\, . . . , /3 n and 7. By Lemma |2.1| for i>2we can write 

^ on = A~ 1 (ei - Si) = af - a~ 

where af and a^~ are in N n . Notice that af 7^ and a~ 7^ 0. Indeed the sum of the entries of 
A^ 1 (ei) is equal to 1/d. Thus \ai\ = \af\ — \a~\ = 0, and consequently the positive and negative 
part of a{ are both non zero for i > 2. The vector af can be written as af = (af, . . . , af n ) for 
i>2. For 1 < k < n consider the integers given by 

m k = max {af, } 

2<i<n lK 

and set /?i = (mi, . . . , m n ). Since f3\ > af, for each i > 2 there is 6i € N n such that f3\ = 6i+af . 
Therefore 

a-i = A~ 1 (ei - ej) = af - a^ = Pi - (9i + a~). 

We set Pi = 9i + a~ for i > 2. Since we have AP\ — e\ = APi — for i > 2, it follows readily that 
APi — e\ > (make i = 2 in the last equality and compare entries). Thus, setting 7 := AP\ — e\, 
it follows that P\,...,p n and 7 satisfy (a) . If each row of B has some zero entry the proof 
of the existence is complete. If every entry of a row of B is positive we subtract the vector 
1 = (1, . . . , 1) from that row and change 7 accordingly so that (a) is still satisfied. Applying 
this argument repeatedly we get a sequence P\ , . . . , p n , 7 satisfying (a) and (b) . 

To see the last part of the assertion, note that condition (a) is equivalent to the equality 
AB = T + 1, where T is the matrix all of whose columns are equal to 7. Since APi = 7 + 
for all i and det(S) = ±det(r + I)/d it suffices to show that det(r + J) = I7I + 1. The latter 
is a classical calculation than can be performed in various ways. One manner is to express this 
determinant as the sum of the determinants of two simpler matrices of the same size and proceed 
by recurrence. One has a more general statement which is given below for lack of an appropriate 
reference with a clear proof: 

Lemma 2.3. Let T = (7$ »•) be an nx n square matrix over an arbitrary commutative ring, and 
let D = diag(c?i, . . . , d n ) be a diagonal matrix over the same ring. Then 

det(r + £>) = det(r) 

+ ^djA[ n ]\{ i} + ^ rf H^2 A [n]\{u,i 2 } 

+ • ' ' + ^2 d h'-- rf i B -i A [n]\{*i,... ) i B _i} 

l<ii<-<i n _i<n 

+ det D, 

where [n] = {1, . . . ,n} and Ayuj^^ ^i denotes the principal (n — k) x (n — k)-minor ofT with 
rows and columns [n] \ . . . , ij.}. 

In particular, ifT has rank at most one and D is the identity matrix, one gets det(r + /) = 
trace(r) + 1, as required above. 



For the proof of the lemma, one notes that, due to the multilinearity of the determinant, the 
following equality of determinants holds: 



■1 
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/7l,l + d\ 71,2 
72,1 72,2 + d 2 



dot 



71, n \ 

72, n 



+ dj 



del 



/7l,l 7l,2 

72,1 72,2 + d 2 



\7»,1 



71, n \ 

72, ™ 



+ d n J 



+ dct 



/rfl 71,2 
72.2 + d 2 



\0 



71, n \ 

72, n 



The recurrence for the second determinant on the right-hand side is obvious, while for the first 
determinant we repeat the procedure, thus writing 



del 



/7l,l 

72,1 
73,1 



71,2 
72,2 + d 2 
73,2 



71,3 
72,3 
73,3 + £*3 



71, n 

72, n 

73,71 



= det 



+ d n J 



+ det 



/7l,l 7l,2 71,3 

72,1 72,2 72,3 
73,1 73,2 73,3 + d$ 



/7l,l 

72,1 
73,1 



71,™ 

72,11 
73.li 



7,1 


2 7n,3 


7»,i 





71,3 


71, 


do 


72,3 


72, 





73,3 + ^3 


73, 



+ dj 
' \ 

+ dj 



□ 



W,l 7„, 3 

and so forth. 

This gives the proof of the lemma and completes the proof of Theorem 12.21 

Remark 2.4. Bridges and Ryser [2] (cf. {4J Theorem 4.4]) studied the equation AB = T + I, 
when A,B are {0, 1} matrices and V is the matrix with all its entries equal to 1. They show 
that if equality occurs, then each row and column of A has the same number r of ones, each 
row and column of B has the same number s of ones with rs = n + 1, and AB = BA. For the 
applications we have in mind here (see Section [3]) this case is uninteresting; furthermore, the 
commutation AB = BA will seldom be fulfilled. 



The proof of Theorem 12.21 provides an algorithm to compute the vectors (3\ , . . . , (3 n and 7 (see 
Example 12.51 and the proof of Proposition 13.51 for specific illustrations of the algorithm) . Other 
means to compute these vectors will be discussed in Section [H 

Example 2.5. Consider the following matrix A and its inverse: 

/d d-l \ 
,4= 1 d — II; A' 
\0 1 



-1 



d 





To compute the /3j's and 7 we follow the proof of Theorem I2.2L Then f3\ = (2,d, 0), 
+ 1,0), p 3 = (d, 1,1), 7 = (d 2 + d- 1, d, 0), and 

/2 1 d> 
B=\d d + l 1 
\0 1, 

By subtracting the vector (1,1,1) from rows 1 and 2, we get 

1 d-l N 
B' = I d - 1 d 
1 
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The column vectors # = (l,d- 1,0), (3' 2 = (0,d,0), $ = (d - 1,0,1), i = (d 2 -d,d- 1,0) 
satisfy (a) and (b). 

3. Application to Cremona maps 

Let R = k[x±, . . . , x n ] be a polynomial ring over a field k and set x a := x^ 1 ■ ■ ■ x® n for 
a = (ai, ■ ■ ■ ,a n ) £ N ra . In the sequel we consider a finite set of distinct monomials F = 
{x Vl , . . . , x Vq } C R of the same degree d > 1 and having no non-trivial common factor. We 
also assume throughout that every Xi divides at least one member of F, a harmless condition. 
The set F defines a rational (monomial) map p n_1 — -> p«-! which will also be denoted F and 
written as a tuple F = (x vi , . . . ,x Vq ). This map is said to be a Cremona map (or a Cremona 
transformation) if it admits an inverse rational map with source P ra_1 . Note that a rational 
monomial map is defined everywhere if and only if the defining monomials are pure powers of 
the variables, in which case it is a Cremona map if and only if d = 1 (the identity map). 

Finally, the integer d is often called degree of F (not to be confused with its degree as a map). 

Theorem 3.1. Let F : P" _1 p™- 1 stand for a rational map defined by monomials of fixed 
degree. If F is a Cremona map then its inverse is also defined by monomials of fixed degree. 
Moreover, the degree as well as a set of monomials defining the inverse can be obtained explicitly 
in terms of the given set of monomials defining F. 

Proof. Let F be a Cremona map defined by a set of monomials fi, ■ ■ ■ , f n of the same degree 
d > 1. By [£l Theorem 1.2] the matrix of exponents of these monomials (i.e., their "log- matrix" ) 
A has determinant ±d. Therefore Theorem 12.21 implies the existence of an n x n matrix B such 
that AB = T + I, where V is a matrix with repeated column 7 throughout. Let gi, . . . ,g n denote 
the monomials whose log-matrix is B and call G the corresponding rational monomial map. 
Letting x 7 denote the monomial whose exponents are the coordinates of 7, the above matrix 
equality translates into the equality 

(h(gi, ■ ■ -,9n), ■ ■ -,fn(gi, ■ ■ -,g n )) = {x 1 ■ x u . . . ,x 7 • x n ). 

Thus the left hand side is proportional to the vector (x±, . . . , x n ) which means that the composite 
map F o G is the identity map wherever the two are defined (see [HI proof of Proposition 2.1]). 
On the other hand, since B is the log-matrix of g\, . . . ,g n , Theorem 12.21 applied in the opposite 
direction says that G is also a Cremona map. Therefore G has to be the inverse of F, as required. 
Finally, notice that the proof of Theorem 12.21 provides an algorithm to compute B and 7. The 
input for this algorithm is the log-matrix A of fx, . . . , f n . □ 

We will call a Cremona map as above a Cremona monomial map. The theorem allows us to 
introduce the following group. 

Definition 3.2. The Cremona monomial group of order n — 1 is the subgroup of the Cremona 
group of P n_1 whose elements are Cremona monomial maps. 

Here we will not distinguish between rational maps F, G : P™- 1 — -> P™- 1 defined, respectively, 
by forms fx , . . . , f m of the same degree and their multiples gx = gfx , ■ ■ ■ -,9m = gfm by a fixed 
form g of arbitrary degree. 

There is a potential confusion in this terminology. For instance, in [6] one allows for more 
general maps by considering the free product with certain Klein groups. Since our goal is 
combinatorial we will stick to the above definition of the Cremona monomial group, but will 
allow for source and target permutation of variables. 
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As a matter of notation, the composite of two Cremona maps F, G will be indicated by FG 
(first G, then F). Likewise, the power composite FF ■ ■ ■ F of m factors will be denoted F m . 

We shift our attention to plane Cremona monomial maps, i.e., we will study the structure 
of the Cremona monomial group for n = 3. Consider the maps H = (x\,x\X2,X2X 3 ) and 

S = (XIX2,XIX 3 ,X 2 X 3 ). 

• For any d > 1, one has H d ~ l = (x d ,x d ~ 1 x 2 ,x d T 1 x 3 ). 

This is a straightforward composite calculation: by induction, we assume that H d ~ 2 is defined 
by x d ~ 1 ,x d - 2 x 2 ,x d - 2 x 3 , hence H^ 1 = H d ~ 2 H is defined by xf dr ~ 1 \ x\ d ~ z x 2 , xf^x^xs which 
is the same rational map as the one defined by x d , x d ~ l x 2 , x d ~ l x 3 (by canceling the gcd x d ~ 2 ) 

• For any d > 1, up to permutation of the variables both in the source and the target, H d ~ l 
is involutive, i.e., coincides with its own inverse - in [10] this was called " p-involutive" . 

Namely, consider the Cremona map G = (xix d T 1 ,x d ,x d ~ 1 x 3 ). Again, one finds straightfor- 
wardly that H d_1 G is defined by xfx^ 1 \x d ~ 1 x 2 d ^ +d x d , x d ~ 1 x 2 i< " d ^x 3 ; canceling the gcd 
xf~ 1 x d ^ d ^ yields the identity map (x\,x 2 ,x 3 ). 

• The monomial maps S and H "nearly" commute; actually they are conjugate by a trans- 
position. 

This is once more straightforward: the general composites SH^ 1 and H d ~ 1 S are defined, 
respectively, by and 3 and these are the same up to 

transposing x\ and x 3 and the extreme terms. This means that they are conjugate by a trans- 
position. It has as a consequence that the subgroup of the Cremona group of P 2 generated by 
S and H is Abelian up to free product with suitable permutation subgroups. 

Lemma 3.3. Let F be a plane Cremona map of degree d of the form F 

Then, up to permutation of the source (variables) and the target (monomials) , F is one of the 
following two kinds: 

F = (xf,x 2 x d ~ 1 ,x d ~ 1 x 3 ) or F = (xix 2 ,x 2 x 3 ,xix 3 ). 

Proof. By hypothesis d = at + a2 = &2 + b 3 = c\ + c 3 = a±b2C 3 + 010263. Then one has 
(3.1) ai(b 2 c 3 - 1) = a 2 (l - c%b 3 ). 

The cases below can be readily verified using these equations. 

Case (I): a% > 1, a 2 = 0. Then a\ = d, 62 = c 3 = 1, and F = (x d ,x 2 x d ~ 1 ,x d " 1 x 3 ). 

Case (II): a% = 0, 02 > 1. Then 02 = d, b 3 = c\ = 1, and F = (x d , x d ~ 1 x 3 , xixf^ 1 ). 

Case (III) (a): 01 > 1, a 2 > 1, b 2 c 3 = 0,b 2 = 0. Then F = (x d ~ l x 2 , x d , xixf' 1 ). 

Case (IH)(b): a x > 1, a 2 > 1, b 2 c 3 = 0, c 3 = 0. Then F = (x d ~ 1 x 2 ,x^ 1 x 3 ,xf). 

Case (III)(c): a 1 > 1, a 2 > 1, c x b 3 = 0, a = 0. Then F = (xix^ 1 , x 2 x d ~ 1 , x d ). 

Case (III) (d): a x > 1, a 2 > 1, c x b 3 = 0, 63 = 0. Then F = (x^' 1 , x%, x(^x 3 ). 

Case (IH)(e): ai > 1, a 2 > 1, b 2 c 3 - 1 > 0, 1 - cib 3 < 0. From Eq. (pTlj) we get &2C3 = 1 and 
c\b 3 = 1. Then F = (x\X2,X2X 3 ,x\x 3 ). □ 

We next give a purely integer matrix theoretic proof of the following result (see [6] for similar 
dealings). 

Proposition 3.4. If n = 3, then up to permutation of the variables (both in the source and 
the target), the plane Cremona monomial group is generated by the maps S,H, where S is the 
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quadratic map of first kind defined by xxx 2 ,xxx 3 ,x 2 x 3 (Steiner involution) and H is the quadratic 
map of second kind defined by x\,xxx 2 ,x 2 x% ( "hyperbolism") . 

Proof. Let F = (x^x^Xg 3 , x^x^x^, x^x^Xg 3 ) be a plane Cremona map of degree d. We 
will prove that, up to a permutation of the variables and the monomials, every plane Cremona 
monomial map is of the form F ■■■ SH d ^ SH dl w ■••(?, where F, G are in {S, H d | d > 1}. The 
proof is by induction on the degree. Consider the corresponding log-matrix 

(ax h cA 
a 2 h c 2 . 

Up to permutation of variables and monomials, there are essentially two cases to consider: 

fax bx 0\ fax cA 

A = \a 2 b 2 or A = I a 2 b 2 . 

\o 3 dj \0 b 3 C 3 J 

By Lemma [3. 3 1 we may assume that A is the matrix on the left, i.e., F = (x^x^Xg 3 , x^x^ 2 , xf). 

Case (I): ax >bx- From ax + a 2 + 03 = 61 + b 2 = d, we get that b 2 > a 2 . Consider the Steiner 

involution given by 5 = (X1X2, X2X3, X1X3). Then 

qp _ (a 1 +b 1 a 2 +b 2 a 3 bx b 2 d ax a 2 _d+a 3 \ _ 61 a 2 a 3 , -ax 62 b 2 -a 2 d-a 3 ax-bx d\ 
lj ± — VI 2 3 1 \ 2 3 ' 1 2 3 / — «a<^ ^2 3 Vl 2 ' 2 3 3 / * 

Thus by Lemma 13.31 we get that ST = F a i+ fe 2-i for some hyperbolism or SF is a Steiner 
involution. Thus multiplying by the inverse of S, we obtain that F has the required form. 

Case (II) (a): bx > ax and det(A) = d. By Lemma 13.31 we may assume ax > 1- From 
the equality d = det(yl) = d(axb 2 — a 2 bx), we get that b 2 > a 2 . Consider the hyperbolism 
H = (x\ , X1X2, £2X3). In this case we have 

fjp /_2ai 2o2_2os „ai+bl x a2+b2„a3 J>x b 2 d\ 

11 1 — v 1 2 3 ' X 2 3^12 3/ 

_ „ai+l„a2„03 f„ai-l 02.03 3>\-X b 2 b\-a\-\ b 2 -a 2 d-a?.\ 

— A± X 2 ^3 I^Xj X 2 Xg jX-^ X 2 X 2 Xg j 

— ™.<H+l™a2„a3 Z7> 

Since iq has degree at most d — 1, we have lowered the degree of i* 1 . Thus by induction F has 
the required form. 

Case (II) (b): bx > ax and det(^4) = — d. We may assume that a 2 > b 2 , otherwise we may 
proceed as in Case (II)(a). By Lemma 13.31 we may also assume that ax > 1. We claim that F 
has the form F = (xf~ 2 x 2 x 3 , xf~ 1 x 2 , x d ). The condition on the determinant is equivalent to the 
equality 0261 — axb 2 = 1. Thus a 2 > 1. Hence using that bx > ax + 1 one has 

a 2 bi > a 2 ax + a 2 > b 2 a\ + 1 = a 2 bx- 

Consequently a 2 = b 2 = 1. The condition det(A) = —d, becomes ax = bx — 1 = d — 2. Therefore 
F has the asserted form. Consider the hyperbolism H = {x\,xxx 2 ,x 2 x 3 ). It is easy to see that 
HF is equal to x 7 (x^ _3 X2X3, x^~ 2 X2, Xg -1 ) for some monomial x 7 . Thus we have lowered the 
degree of F and we may apply induction. □ 

The next result is classically well-known. One has a simple direct proof in the case of Cremona 
monomial maps. The proof shows moreover explicit simple formulae for the Cremona inverse in 
the case of 3 variables. 

Proposition 3.5. If n = 3, then a Cremona monomial map and its inverse have the same 
degree. 




8 ARON SIMIS AND RAFAEL H. VILLARREAL 

Proof. The proof is based on the method employed in the proof of Theorem 12.21 Let A be the 
matrix defining a Cremona map of degree d: 

/ ax bi ci 
\a 3 h c 3 , 

For simplicity we assume that det(A) = d, the case det(A) = —d can be shown similarly. Up to 
permutation of variables, as in the proof of Proposition 13.41 there are essentially two cases to 
consider: 

( ai b x 0\ 
A = \a 2 b 2 or A 

\a 3 d) 

It is readily seen that the inverse of A is given by 

/ b 2 -bx \ 

A' 1 = -a 2 ax or A' 

y-a 3 6 2 /d bxa 3 /d 1/dJ 

respectively. By the argument given in the proof of Theorem 12.21 we get that 

0x = (d,O,Q), 7 = (dax - l,da 2 ,da 3 ) or 0x = (b 2 ,0,b 3 ), 7 = (axfa + b 3 cx - 1, a 2 b 2 , b 3 c 3 ) 

respectively. Therefore 

id bx \ (b 2 cx 0' 

B = ax + a 2 a 2 or B = c 3 a 2 

\0 a 3 (a 3 b 2 + l)/d) \b 3 a X/ 

respectively. Using that the sum of the entries in every column of A is equal to d and det(vl) = d, 
it is seen that the sum of the entries in each column of B is equal to d and det(i?) = d. □ 




4. Further combinatorial and computational aspects 



Let vx, ■ ■ ■ , v n be a set of vectors in N n such that \vi\ = d > 1 for all i and det(^4) = ±d, where 
A is the n x n matrix with column vectors v\, . . . , v n . Then, by Theorem 12.21 there are unique 
vectors . . . , (3 n , 7 € N ra such that the following two conditions hold: (a) Afc = 7 + for all 
i, and (b) The matrix B whose columns are . . . , j3 n has at least one zero entry in every row. 

The proof of Theorem 12.21 provides an algorithm to compute 0x, . . . ,/3 n , 7. Next we discuss 
two other methods, one of them based on integer programming techniques and the notion of a 
Hilbert basis. To compute the sequence /3\, . . . , f3 n , 7 using linear programming we regard the 
and 7 as vectors of indeterminates and introduce a new variable r. Consider the homogeneous 
system of linear inequalities 



(4.1) Afc = 7 + rei, i = l,...,n 

Pi > 0, i = 1, . . . ,n 
7 > 0, r > 0. 

This linear system has n 2 equality constraints and I = n 2 + n + 1 indeterminates. The set C 
of solutions form a rational pointed polyhedral cone. By Theorem 16.4], there is a unique 
minimal integral Hilbert basis 

H = {hx, ■ ■ ■ , h r } 
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of C such that Z^nR+7Y = NTL and C = W + TC (minimal relative to taking subsets), where W + Ti 
denotes the cone generated by TC consisting of all linear combinations of TC with non-negative real 
coefficients and N7Y denotes the semigroup generated by TC consisting of all linear combinations 
of TC with coefficients in N. The Hilbert basis of C has the following useful description. 

Theorem 4.1. [7J p. 233] TC is the set of all integral vectors ^ h € C such that h is not the 
sum of two other non-zero integral vectors in C. 

Proposition 4.2. There is a unique element h of TC with r = 1 and this element gives the 
unique sequence Pi, . . . , f3 n , 7 that satisfies (a) and (b). 

Proof. We set xq = (Pi, . . . , (3 n , 7, 1) £ j^™ 2 +™+ 1 5 where Pi, . . . , /3 n , 7 is the unique sequence that 
satisfies (a) and (b). First we show that xq is in TL. Clearly xq is in C. Thus we may assume 
that xq is written as 

xq = m h i H r- rj k h k , ^ r/j G N for i = 1, . . . , fc, 

where /i^ has its last entry equal to 1 and the last entry of hi is equal to for i < k. The vector 
hk has the form 

where the 0^ s and 7^ are in N n and satisfy Eq. (14. ip . i.e., they satisfy (a). Notice that from 
the first equality one has xq > h^. Then Pi > for all i and 7 > 7^). Therefore the P^'s 
and 7" also satisfy (b). Consequently xq = h^ and xo S TC, as claimed. Let hi be any other 
element in Ti whose last entry is equal to 1. Next we show that hi must be equal to xo- The 
vector hi has the form 

h t = (pf\...,P^,^\l), 

where the Pj^'s and 7® are in N ra and satisfy Eq. (|4.1h . Since the Pj's and 7^ satisfy (a), it 
suffices to show that they also satisfy (b). We proceed by contradiction. Assume that the Ith 

(i) 

entry of P - is not zero for j = 1, . . . , n. For simplicity of notation assume that i = 1. Then the 
vectors 

ti = (ei, . . . , ex, Ae x , 0) and ft" = - e 1; . . . - e 1)7 « - Aei, 1) 
are integral vectors that satisfy Eq. (|4.ip . i.e., ft' and /i" are integral vectors in C and /i = h' + h", 
a contradiction to Theorem 14.11 □ 

There are computer programs that can be used to find integral Hilbert basis of polyhedral 
cones defined by linear systems of the form x > ; A'x = 0, where A' is an integral matrix. We 
have used [3] to compute some specific examples with this procedure: 

Example 4.3. Consider the following matrix 

/2 1 N 
,4= 1 1 

V° h 

To compute Pi , P2 , P3 , 7 we use the following input file for Normaliz [3] 
9 

13 

210 000 000 -100 -1 
011 000 000 0-10 
001 000 000 00-1 
000 210 000 -100 
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Part of the output file produced by Normaliz is: 

4 generators of integral closure: 
0100100101100 
0010010010110 
1001001002000 
1100201012101 

From the last row we get: X = (l,l,0),/3 2 = (0,2,0),/% = (1,0,1), 7 = (2,1,0). 

Another method to compute Pi, . . . ,(3 n is based on the following result whose proof we omit. 
Of course, to use this result effectively one must seek methods to compute 8. 

Proposition 4.4. Let vx,...,v n be a set of vectors in N n such that \vi\ = d > 1 for all i 
and det(vl) = ±d, where A is the n x n matrix with column vectors v\, . . . ,v n . Then there is 
0^6= (Si) G Ql such that 5 + A~ 1 e l G N n for all i. If 

Pi = A~ 1 (5ivi H ^5 n v n + ei), i = l,...,n 

and 7 = 5±vi + ■ • • + S n v n . Then . . . , /3 n , 7 are vectors in N n that satisfy condition (a). 
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